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We study phase transition of a nonequilibrium statistical-mechanical model, in which two degrees 
of freedom with different time scales separated from each other touch to their own heat bath. A 
general condition for finding anomalous negative latent heat recently discovered is derived a from 
thermodynamic argument. As a specific example, phase diagram of a spin- lattice gas model is 
studied based on a mean-field analysis with replica method. While configurational variables are 
spin and particle in this model, it is found that the negative latent heat appears in a parameter 
region of the model, irrespective of the order of their time scale. Qualitative differences in the phase 
diagram are also discussed. 
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I. INTRODUCTION 

Phase transitions under non-equilibrium conditions 
have attracted with a great deal of attention in 
statistical- mechanical problems [l|, [2j 3 • There have been 
many investigations on non-equilibrium phase transitions 
so far [J, la, la, III I which have revealed a fascinating new 
transition behavior different from equilibrium transition. 
In general, probability distribution in non-equilibrium 
cannot be expressed in terms of only energy functional, 
which causes a difficulty in theoretical study. 

Recently another class of non-equilibrium systems that 
exhibits a phase transition has been studied[8|, in which 
two different degrees of freedom coupled to their own heat 
baths interact with each other through multi-body inter- 
actions. For simplicity, time scales of these two variables 
arc assumed to be well separated. Then, the systems con- 
sist of slow and fast variables belonging to different time 
hierarchy. The fast variables behave in quasi-equilibrium 
for a given set of slow variables that plays a role as 
quenched variables for the fast ones. Meanwhile, the 
slow variables are not given by independent distribution 
function as in quenched disordered systems but affected 
through mean force of fluctuating fast variables. Such 
systems with a hierarchy in scpcrated time scales and 
different heat baths arc here called "two-temperature" 
system. These systems are adopted to describe neural 
network s^tems with a synaptic evolution[9|, evolving 
networks |10l|. and some kind of NMR systems [ll|. In 
contrast to most of non-equilibrium systems, the steady- 
state distribution of the models is formally expressed in 
terms of the energy function using the replica method 
that is a standard tool for studying thermodynamic prop- 
erties of the quenched disordered systems |12|. 

The replica formalism for two-temperature systems has 
been introduced in Refs. l3, ll^, llJl- While the quenched 
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disordered systems require the replica limit in which the 
replica number takes to be zero, the two-temperature sys- 
tems have a physical meaning for any value of the replica 
number which corresponds to a ratio of two tempera- 
tures. In this sense, this is regarded as a generalization 
of quenched systems and is sometimes called as partial 
annealing systemQ. A most crutial difference from the 
quenched systems is that the slow variables are also dy- 
namical coupled to their heat bath. Therefore, both the 
slow and fast variables are responsible for phase transi- 
tion. This could provide a new cooperative phenomenon 
over different time scale. 

In fact, AUahverdyan and Petrosyan, hereafter referred 
to as AP, studied a mean-field spin model as a two- 
temperature system and found that the model exhibited 
a first-order phase transition with anomalous negative 
latent heat, which never occurs in equilibrium statisti- 
cal mechanics. However, this peculiar behavior observed 
in the two-temperature system is not well understood. 
We pursue phase transition in the two-temperature sys- 
tems and give a general condition that the system ex- 
hibits the negative latent heat with the help of the idea 
of thermodynamics. It is also found that, in the sys- 
tems, two different entropys associated with the fast and 
slow variables respectively play a competitive role in dc- 
terming the phase boundary of first-order transition. We 
further studied two-tempraturc version of a spin-latticc- 
gas model, similar to that studied by AP, as a specific 
example. The spin-lattice-gas model consists of two de- 
grees of freedom, spins and particles, which has been 
studied for a given Hamiltonian in equilibrium ll5|. The 
two-tcmpcraturc version is characterized by not only the 
Hamiltonian but also the order of time scales of two vari- 
ables. AP studied the case where the spins were slow and 
the particles were fast. We study this case with some 
modified Hamiltonian and also the other case, namely 
the spins and particles behave as fast and slow variables, 
respectively. We then find that the existence of the neg- 
ative latent heat is common to both cases, suggesting 
that it is observed in a wide class of the two-temperature 
systems. On the other hand, qualitatively different be- 



havior is also found in phase diagram of the two cases, in 
particluar, the stabihty of ferromagnetic ordered phase. 
This paper is organized as foUows. In Sec. |lll we re- 
view the rephca formaUsm of the two-temperature sys- 
tem, which leads to an equilibrium model of a replicated 
system. We also discuss phase boundary of first-order 
transition and derive a Clausius-Clapeyron relation in 
this system. This relation enables us to find generally a 
geometric property of the phase boundary and the neg- 
ative latent heat. In Sec. Illli we explicitly define two 
mean-field spin-lattice-gas models with different order of 
time scales and give self-consistent equations for the mod- 
els. The results obtained by solving the equations are 
presented in Sec. IIVI In Sec. |Vl we summarize our re- 
sults. 



II. TWO-TEMPERATURE FORMALISM WITH 
DIFFERENT TIME SCALES 

In this section, we review a theoretical formalism for 
two-time-scale and two-temperature system[8|. Suppose 
a system described by a Hamiltonian H{f,s), in which 
/ is a symbolic notation of fast degree of freedom and 
s is of slow degree of freedom. These variables / and s 
are in contact with their different heat bath with tem- 
perature Tf and Tg, respectively. We assume that the 
two characteristic time scales on the variables s and / 
are well separated from each other and that the thermal 
average of the fast variable / can be taken with a fixed 
configuration of the slow variable s. Then, the condi- 
tional probability P(/|s) of finding a configuration / for 
a given s at the inverse temperature /3/ = l/Tf is defined 
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where the normalization constant or the partition func- 
tion of the fast variable is set as 



Z(s) = TYfe-^f"^f^'^ 



(2) 



Hereafter the Boltzmann constant is set to be unity. One 
can define partial free energy for the fast variable as 
Ffis) = -Tf\ogZis). 

The steady-state probability of slow variables P{s) 
is derived by an adiabatic approximation of two- 
temperature Langcvin equation [l3[ . The force acting on 
s is assumed to be an averaged derivative of the Hamilto- 
nian with respect to the slow variable over the conditional 
probability, which is represented by the partial free en- 
ergy as d^- The equilibrium distribution P{s) at 

the inverse temperature Ps = \/Tg is given by 



P{s) = 



.-f>.Fi[s) 



where 



T.e-^^^^(^). 



(3) 



(4) 



The total free energy T is defined by .F = —Tg log-E. 

Using the replica trick, the model can be mapped onto 
an equilibrium problem with a replicated Hamiltonian for 
the integer number of ratio n = Tf /Tg, 
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where J^'-* denotes replicated fast variables. This could 
be extended to any real value of the ratio Tf/Tg after 
calculating the replicated system in a standard man- 
ner of the replica method. While one takes the replica 
limit 77. ^ for the quenched disordered system like 
spin glasses, any value of n makes a sense as the two- 
temperature system in this context. This formalism is 
also interpreted as a kind of statistical-mechanical prob- 
lems with randomness. In particular, note that the dis- 
tribution of the random variables is determined by not 
only a given independent function but also the thermal 
averaged quantity of the fast variables. The latter leads 
to a non-trivial correlation among the slow variables. 

We discuss thermodynamic properties of the two- 
temperature system. The simultaneous probability 
PU,s) is expressed as P(/, s) = P{f\s)P{s). The to- 
tal entropy S defined by the simultaneous probability is 
decomposed into two degrees of freedom as 

S = -TisjPif, s) logP(/, s)=S, + Sf, (6) 

where Ss and Sf are expressed as 

S, = -Tr,P(s)logP(s), (7) 

Sf = -TT,P{s){TTfP{f\s)logP{f\s)). (8) 

The total free energy J- is formally expressed as 



TiTs, Tf) = {H{f, s))f - TfSf - T^Ss, (9) 

where (■■•)/ and '^^ denote an average over the variables 
/ with P{f\s) and s with P{s), respectively. It should 
be noted that the free energy is also rewritten by 



HTs,Tf)^Ff{s)-T,Ss. 



(10) 



Averaging over the fast variables /, the thermodynamic 
structu re is fo und by regarding the averaged partial free 
energy Ff{s) as an "energy" for the slow variable s. 
Namely, the averaged partial free energy and the entropy 
Sf for the slow variables decrease monotonically with de- 
creasing Tg. 

As a consequence of the thermodynamic structure [l6|, 
a Clausius-Clapeyron like relation for two-temperature 
systems is derived, which gives us a topological property 
of a first-order-transition line. Suppose a phase diagram 
of the system onto two-temperature plane of Tg and Tf. 
We take two points {Tf, T,) and {Tf+STf,Ts+STs) which 
are located on either side of the first-order-transition line. 
The free-energy difference SF between these points with 
small temperature differences STf and STs is given by 



SJ^-^ 



SfSTf - SgSTs 
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At the first-order transition point (Ti ,Ti '), the or- 
dered and disordered states coexist and the free energy 
of these states coincides with each other, meaning 

AT = J^(°)(T^^', Tji)) - T'-''\T'f^\TP) ^ 0, (12) 



where the upper sufSxes o and d denote the ordered 
and the disordered states, respectively, and A A means 
the difference of a physical quantity A between the or- 
dered and disordered states at the transition point. Us- 
ing Eqs. (fTTj) and (fT^ . the Clausius-Clapevronflq like 
relation is obtained as 
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This implies that when the slope of the phase boudary 
dTs /dTr is positive, ASs and ASf arc opposite to 
each other. The frcc-cnergy difference AT is also ex- 
pressed as AT = AW - TfASf - T^ A5s with U being the 
internal energy {H)f. Thus, we obtain the relation be- 
tween the deference of the internal energy and the phase 
boundary as 



AU{Ty\TP) = tPASs il-n 




(14) 



Because the entropy Sg is a monotonically decreasing 
function of Ts , the sign of AU depends on only the gradi- 
ent of the phase boundary. This implies that the condi- 
tion to find the negative latent heat is ^ < dTg /dTr 
when Ts decreases. On the other hand, when Tf de- 
creases, the condition for the negative latent heat changes 
to < dTP/dTJ:^^ < i. While AP exphcitly found that 
a specific spin-lattice gas model exhibited the negative 
latent heat in a region of the phase diagram using the 
replica method, wc find a general condition for which 
the negative latent heat appears thorough the thermo- 
dynamic argument. 



III. MEAN-FIELD SPIN-LATTICE GAS MODEL 

A model Hamiltonian we studied is an infinite-range 
spin-lattice-gas model, that is given by 

, N 

Hi{S^,p.}) = --J2(JS^S,+e^p,p,+aJ2p^ (15) 



m 



i=l 



where 5*^ = ±1 arc spin variables, pi — 0,1 are parti- 
cle occupation variables and they are defined on N sites. 
In the case where the spins 5*^ are the slow variable and 
the particles pi the fast, referred to as case-1 model, the 
model system with ej = is identical with that studied 
by AP[^. We also consider the inverse case where the 
spins and the particles represent the fast and slow vari- 
ables, respectively, which is referred to as case-2 model. 



The spin and particle variables are coupled to their own 
heat baths whose temperature is denoted by Ts and Tp, 
respectively. The sum is taken over all pairs of sites. The 
interactions J and ef denote a ferromagnetic coupling 
and an attractive interaction between particles, respec- 
tively. In this paper, J is taken as a unit of energy and 
temperature. The first term of the Hamiltonian consists 
of a spin mediated interaction and direct one. The sec- 
ond term plays a role for controlling a particle number 
with chemical potential a which is chosen to be a posi- 
tive value. The spin-lattice-gas model could exhibit two 
types of phase transition which are magnetic and den- 
sity orderings. The interaction —yJSiSj + e/ ) between 

particles prefers to increase the particle density and mag- 
netically ferromagnetic ordering, while the chemical po- 
tential a tends to decrease the particle density. In this 
sense, these two energy terms compete with each other. 
Furthermore, two different kinds of the entropy associ- 
ated with the fast and slow variables also compete with 
the energy terms. 

Since the Hamiltonian is an infinite range model, the 
trace of Eq. ([5]) is carried out with the help of the replica 
method by introducing two auxiliary fields m and p, 
which correspond to average magnetization and particle 
density, respectively. The self-consistent equations for m 
and p are written as, for the case-1 model. 



^ S(j){S;m,p){l + (j){S;m,p)Y' 



s=±i 



Y, (l + 0(S;m,p))" 

s=±i 

^ (j){S\m,p){l + (j){S\m,p))'"'~ 
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s=±i 



Y, (l + 0(5;m,p))" 



(17) 
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and for the case-2 model, 
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2"+( E ^iS;m,p)) 
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where 0(5"; m, p) — e V / . Here /Sf is the 

inverse of the fast temperature, which corresponds to 
1/Tp in the case-1 model and l/Ts in the case-2 model. 
The free energy of system is represented with a solution 
(?7io, po) of the above self-consistent equations as 



J^i{mo,po,Ts,Tp) 



JniQ 



e/Po 



0.2 



S 0.1 



0.05 



0.072 


' 


' 




' ' 


\ ' 


' 


' 


EpO.O 

y 

y 
y 
y 

y 








r 


S 0.071 






/ 


y 

y 


PM 




V 


U 


y 
y 
y 
y 
y 
y 


0.07 










0.1 


0.15 


0.2 


y 






I'P 


/ 
/ 
/ 
/ 






FM 











0.05 0.1 0.15 0.2 0.25 0.; 

1p 



art 



FIG. 1: Phase diagram of the case-2 model on Ts — Tp plane 
with a — 0.45 and e/ = 0. PM and FM denote paramagnetic 
and ferromagnetic phases, respectively. The solid and dashed 
lines represent the first and second order transitions and the 
dotted one represents the instability limit of the paramagnetic 
solution. The inset shows an enlarged view around the critical 
point. 



Tsiog( Y, {i + HS; 

\S=±1 



mo,pa)y 
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for the case 1 and 



J^2{mo, po, Ts, Tp) = - yJm^ + ejplj 



Tplog(2" + ( ^ ^{S; 
\ s=±i 



mo,Po)r] (21) 



for the case-2. Here n is defined as the ratio of the fast 
temperature to the slow one, namely n is tt^ and 7# for 

^ s ^ p 

the case-1 and casc-2, respectively. In this paper, jFi and 
^2 as functions of m and p are loosely called to "free 
energy" in a sense of Ginzburg-Landau free energy. Ac- 
cording to the Eqs. ([7]) and dH]), two kinds of decomposed 
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FIG. 2: Tp dependence of U, TsSs, and averaged partial 
U - TsSs at a = 0.45 and Ts = 0.071 
There are three phase transitions, in- 



free energy Fs{p) = 
in the case-2 model 
dicated by vertical lines, with Tp changing for a fixed Tg. 
Two first-order transition occurs at Tp = 0.116 and 0.180. 
The former is between PM and FM, whereas the latter is a 
re-entrant transition between dense and dilute ferromagnetic 
phases. Between these transitions, a second-order transition 
occurs at Tp — 0.175. 



transition involves the negative latent heat discussed in 

secnn 

We study phase diagram of the case-2 model, the time 
scale reversed version studied by APQ. Figure [T] shows a 
phase diagram on the Ts — Tp plane for the case-2 model 
with e^ = and a = 0.45, which is in comparison to 
the phase diagram of the case-1 model shown in Ref. Q 
under the condition of ey = 0. 

While the first-order phase-transition temperature 
shows rather weak dependence of Tp and takes a finite 
value at Tp = 0, it behaves non-monotonically as a func- 
tion of Tp near the critical point as shown in the inset 



of Fig. [T] According to Eq. p^ , in the region between 
A and B shown in the figure, the latent heat becomes 
anomalously negative when Tp decreases. Fig. [5] shows 
Tp dependence of thermodynamic quantities for a fixed 



entropy are termed as Ss and Sp, respectively for the Ts, where phase transitions occur three times as a func- 
case-1 model, and Sp and Ss for the case-2 model. tion of Tp. As Tp decreases at Ts = 0.071, a first-order 

transition occurs at Tp = 0.18 from a dense ferromag- 



IV. RESULTS AND DISCUSSIONS 

A. Phase diagram and negative latent heat 

We first discuss phase diagram of the spin-lattice gas 
model with the condition £/ = both for the case-1 and 
the case-2 models. The case-1 model with Cf = 0, that 
is the same as that studied by APQ, shows that a ferro- 
magnetic phase has a place in a low Ts region at a = 0.45, 
and that there is a region of the phase boundary in which 
the internal energy of the ferromagnetic phase is higher 
than that of the paramagnetic phase. Namely, the phase 



netic phase to a dilute one and a second-order transition 
between the dilute ferromagnetic and the paramagnetic 
phases at Tp = 0.175. Eventually, the transition from 
the paramagnetic to the ferromagnetic phases again oc- 
curs Tp = 0.116 

At the highest transition temperature Tp ~ 0.180, the 
internal energy and the entropy Ss have a positive jump, 
while the averaged partial free energy Fsip) decreases 
monotonically. This means that the internal energy of 
highly ordered phase at higher temperatures is lower than 
that of disordered phase at lower temperatures. A sim- 
ilar first-order transition is found between A and B in 
Fig. [1] Thus, this phase transition could be simply in- 



tcrpreted as a kind of reentrant transition, in which the 
low-temperature disordered phase is stabihzcd by an en- 
tropic effect. This is in contrast to the casc-1 model 
where the low-temperature phase is the disordered para- 
magnetic one. 

Another difference between the case-f and case-2 mod- 
els is found in the topology of the phase diagram. 
Whereas the case-1 model has a tricritical point at which 
the first and the second-order-transition lines merge, the 
first-order-transition line enters into the ferromagnetic 
phase in the case-2 model as shown in Fig. [TJ Inter- 
estingly, a density origin phase transition occurs in the 
ferromagnetic phase. Near the transition the free-energy 
has four different local minima which correspond to high- 
density and low-density ferromagnetic states and their 
time-reversal ones. This is qualitatively different from 
that observed by AP in the case-1 model. 

Let us discuss the effect of e/ term in the spin-lattice 
gas model. The first-order transition of this system is 
originated with the particle density. Therefore, the first- 
order-transition line could be changed by introducing 
the direct interaction between particles, the ey term in 
Eq. (fTS]) . We study how the effect of the e/ term on 
the phase diagram of both the case-1 and case-2 model. 
First, wc focus on ej^-dcpendcnce of the region in which 
the negative latent heat is observed. Fig. [3] shows the 
phase diagram with ey = 0.4 in the casc-I model and 
the inset shows that with e/ = 0. As the value of e/ in- 
creases from zero, the first-order transition temperature 
Tj for a fixed Tp increases and the ferromagnetic region 

is extended. Non-monotonic behavior of Tg found in 
the inset of Fig. [3] near the multicritical point disappears 
with e/ increasing. Eventually, at the value e/ = 0.40 as 
shown in Fig. [31 the first-order-transition line is mono- 
tonic as a function of Tp. The argument in Sec. [IT] yields 

that the monotonic behavior of Tg as a function of Tp 
means the absence of negative latent heat on the transi- 
tion. Thus, it is found that the region in which negative 
latent heat observed is robust against an infinitesimal 
attractive interaction and disappears by further increas- 
ing the interaction. This suggests that the negative la- 
tent heat is not peculiar behavior in the two-temperature 
system and could be observed by controlling the model 
parameter. Similar behavior is observed in the case-2 
model. The phase diagram with e/ = 0.4 for the case-2 
model is shown in Fig. 21 As seen in the case-1 model, 
the ferromagnetic phase transition is also enhanced and 
the non-monotonic region of the first-order-transition line 
becomes narrow with e/ increasing. 



B. Stability of ferromagnetism in the two models 

In this subsection, we discuss phase diagram with rel- 
atively large ey. There is remarkable difference in the 
stability of ferromagnetic order between the case- 1 and 
case-2 models. Fig. [51 shows a phase diagram of the case- 
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FIG. 3: Phase diagram of the case-1 model on Ts — Tp plane 
with a — 0.45 and e/ — 0.4. The thick line represents the 
first-order phase transition. The symbols of lines are the same 
as those in Fig.l The inset shows the phase diagram of the 
case-1 model with a = 0.45 and e/ = 0. 
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FIG. 4: Phase diagram of the case-2 model with a — 0.45 
and e/ = 0.40. The symbols of lines are the same as those in 
Fig. [l] The inset shows an enlarged view near the multicrit- 
ical point. The first-order-transition line intersects with the 
paramagnetic instability line. 



1 model with ef = 0.8, in which the ferromagnetic phase 
exists stably up to extremely high temperature. As ef 

increases further, Tg takes a finite value in the limit 
Tp = again, as shown in Fig. [51 Namely, in a finite 

range of ef, Tg diverges as a function of Tp and then 
the ferromagnetic phase becomes stable even in the high 
Ts limit. In contrast, the ferromagnetic phase boundary 
in the case-2 model changes modestly with ef increasing 

as shown in Fig. JH and Tg remains to be finite in the 
limit Tp = . This suggests that the difference of the 
time scales between the particle and the spin strongly 
affects on the stability of the ferromagnetic phase. 
It would be helpful for making clear the issue men- 



tioncd above in the two models to see an instability con- 
dition of the paramagnetic phase. The paramagnetic in- 
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on the Ts 



Tp plane is 
-0, 



stability line, (T^' 

simply determined by the condition |^ 

m— 0,p— Pp 

because ofF-diagonal term of a Hessian matrix of the free 
energy with respect to m and p vanishes in the param- 
agnetic phase. Then, the self-consistent equations for p, 
Eqs. (fT7|) and (fT9|) , in the paramagnetic phase are simply 
reduced to the equation. 



(pra) 



p/3p(e/p<,'""'-a) 



1 -|_ e/3p(e/Po " -") 



(22) 



When Ts increases to infinity, Pg " decreases gradually 
down to 1/2, but never reaches to zero. Consequently, 
the magnetization m remains finite even at Tg = oo. 
In fact, in the limit Ts —>■ oo, the free-energy difference 
between the ferromagnetic and the paramagnetic solution 
takes the form — — g-^ + a/2, that is the internal energy 
for the ferromagnetic solution. This yields the stability 
condition of the ferromagnetic phase as ej > 4a — J. 
For example with a — 0.45 and e/ — 0.8, as shown in 
Fig. ini the ferromagnetic phase is extended up to very 
high Ts temperature, although the instability line of the 
paramagnetic solution goes down to the origin. 

For sufficiently large e/, the paramagnetic solution is 
qualitatively changed by the effect of the attractive intcr- 



^(p™) 



where Pq denotes a solution of the equation in the 
paramagnetic phase. By using the solution of the equa- 
tion, the instability condition for the case-1 model is 
given by 



action. Then, p^ 



(pm) 



R Ti ^ o(p'")n ^('""h „(p™)' 



T„ 



Ts' 



= 0. (23) 



-n(pmi) 



This yields the instability temperature Tg as a func- 
tion of Tp expressed as 



1 at Ts — > oo in the limit Tp = 

and the free-energy difference is modified to g-^ — ^ . 

The dense paramagnetic solution becomes dominant at 
{Tp,Ts) = (0, oo). Thus, the first-order transition tem- 
perature Tg (Tp) can diverges only in a finite range of 
ef in the case-1 model. 

In the case-2 model, on the other hand, the spin vari- 
ables fluctuate as fast degree of freedom for a given slow 
particle configuration. The paramagnetic instability con- 
dition is then given by 
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When the denominator 1 
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where Pq ™ is again determined by Eq. 

(pm)^ 



(27) 



The r„ 



IS zero, dependent term coupled to Pq cancels out because of 
j^U^rrn) jj^gj.ggg g^j^j j^j,j^^j, i-jjj, ferromagnetic phase be- the symmetry of the fast spin variable. In the paramag- 
comes stable even at Ts = oo. As a trivial example, netic phase, the particle density pg of the case-2 model 

is the same value as the case-1 model. Thus, T^"^^' could 
not diverge in any value of e/ and Tp, in sharp contrast 
to the case-1 model. This is, however, an necessary con- 
dition but not the sufficient one for the finite transition 
temperature at Tp = 0. 

We see again the phase boundary at Tp = 0. In the 
case-2 model, the only particle configuration that min- 
imizes the partial free energy at Tp = contributes to 
the ensembles and the fast spin variables ffuctuate under 
the resultant particle configurations. The self-consistent 



when Ef = 2a, T^p""'' goes to infinity at Tp = 1/4. At 
(e/,a) = (0.90,0.45), the first-order-transition line and 
the paramagnetic instability line almost merge and the 
jump of thermodynamic quantities at first-order transi- 
tion is quite weak in a wide region of (Ts, Tp) plane. Be- 
cause the instability line is located on the second-order 
transition or inside the ferromagnetic phase, the diver- 
gence of Tg™'^' (Tp) means the stability of ferromagnetic 
phase at an infinite Ts- 

We show explicitly the stability of ferromagnetic phase 
in the case-1 model at Tp = 0. In the case-1 model, the 
spins that arc slow variable can fluctuate even at Tp = 0. 
The particle configuration is adaptively determined for 
a given slow spin configuration by minimizing the free 
energy. For intermediate e/, which is, to be precise, given 
by £/ < 1.45 at a = 0.45, the paramagnetic state is an 
empty state at Tp = 0, namely mg = and ppm = 0. On 
the other hand, the self-consistent equations, Eq. p^ and 

Eg. lfTT)) . for the ferromagnetic solution, nia and Pg " , at 
Tp = are then 



at Tp =: 0, while the 



Tno — ±Pg 



Po 



g/9s((J+e/)p^^"'+a) 



(25) 

(26) 



equation for Pg "" leads to Pg "* 
corresponding equation for the ferromagnetic phase leads 
to a fully occupied solution with pg = 1. For the latter, 
the magnetization tuq is determined by the equation 

niQ = tanh (3s JniQ, (28) 

under the condition e'^^^'^^""^ cosh/?s jTTig > 1. Thus, 
Tg never diverges and the ferromagnetic phase only 

emerges at most Ts < 1/J. Actually, T^ ' is obtained 
by solving the equation 

= ^{Jml + ef) - Ts log {e^s{ef+a) cosh/3s Jmg),(29) 

which is derived from the condition that the free-energy 
difference becomes zero at the transition temperature. 
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FIG. 5: Phase diagram of the case-1 model with a = 0.45 
and €f — 0.80. The symbols of lines are the same as those in 

Fig.m 
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FIG. 6: Phase diagram of the case-1 model a — 0.45 and 
£/ — 0.99. The symbols of lines are the same as those in 

Fig.m 



SUMMARY 



systems, which enables us to link the topology of phase 
diagram and discontinuity of thermodynamic quantities 
at first-order transition. In particularly, a general con- 
dition to find the anomalous negative latent heat that is 
found in a specific spin model[8j] is reduced to a simple 
topological constraint on the phase diagram. To be con- 
crete, when the slope of the first-order phase boundary 
is a certain value determined by the retio of two tem- 
peratures, the negative latent heat appears. It should 
be worth noting that this criteria can be applied to any 
model including short-ranged models in finite dimen- 
sions. 

We have also performed a mean-field analysis of two- 
temperature version of a spin-lattice gas model that has 
spins and particles as configurational variables. Gener- 
ally, two-temperature systems are characterized by the 
Hamiltonian and time-scale order of two varibles. Even 
in the same Hamiltonian, phase diagram still depends 
on the choise of the time-scale order. We have studied 
phase diagram of the spin-lattice gas model for two differ- 
ent cases; one is that the spins arc slow and the particles 
are fast, which is the same as that studied by APQ, and 
the other is alternative. Furthermore, the effect by in- 
troducing preferentially an attractive interaction for one 
of the two variables is studied. We have found that the 
general condition for the negative latent heat is satisfied 
in a parameter region both for two cases, suggesting that 
the negative latent heat is not accidental but frequently 
observed in two-temperature systems. By increasing the 
attractive interaction, the parameter region becomes nar- 
row in common. On the other hand, qualitatively differ- 
ent properties are found in the phase diagram, such as 
the stability of the ferromagnetic order and the existence 
of the ferromagnetic-ferromagnetic transition. This indi- 
cates that the time-scale order plays a significant role in 
phase transitions and cooperative phenomena. An inter- 
esting and open problem would be to see if the results 
found in the spin-lattice gas model are preserved beyond 
the mean-field analysis, for instance in finite-dimensional 
short range models. In this direction, we further progress 
for the model up to the Bethe approximation [17| . 



We have studied phase transition of a non-equilibrium 
statistical-mechanical model that consists of two degrees 
of freedom with different time scales and heat baths, 
called two-temperature systems. A theoretical frame- 
work based on the replica method and its thermody- 
namical structure, which have been already given in the 
literature!^, l^j 13 1 ^^^ summarized. As a direct conse- 
quence of the structure, we have pointed out the existence 
of a Clausius-Clapeyron like relation in two-temperature 
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